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ABSTRACT 

O ! 

, We have cross-correlated X-ray catalogs derived from archival C handra ACIS obser- 

vations with a Sloan Digital Sky Survey (SDSS) Data Release 2 (DR2) galaxy catalog 
to form a sample of 42 serendipitously X-ray detected galaxies over the redshift interval 
0.03 < z < 0.25. This pilot study will help fill in the "redshift gap" between local X-ray- 

00 ' studied samples of normal galaxies and those in the deepest X-ray surveys. Our chief 

purpose is to compare optical spectroscopic diagnostics of activity (both star-formation 

,-h . and accretion) with X-ray properties of galaxies. Our work supports a normalization 

value of the X-ray-star-formation-rate (X-ray-SFR) correlation consistent with the lower 
values published in the literature. The difference is in the allocation of X-ray emission 
to high-mass X-ray binaries relative to other components such as hot gas, low-mass 
■ X-ray binaries, and/or AGN. We are able to quantify a few pitfalls in the use of lower- 

resolution, lower signal-to-noise optical spectroscopy to identify X-ray sources (as has 

"--^ ■ necessarily been employed for many X-ray surveys). Notably, we find a few AGN that 

likely would have been misidentified as non-AGN sources in higher-redshift studies. 
However, we do not find any X-ray hard, highly X-ray-luminous galaxies lacking optical 
spectroscopic diagnostics of AGN activity. Such sources are members of the "X-ray 



CO 



Bright, Optically Normal Galaxy" (XBONG) class of AGN. 



Subject headings: diffuse radiation — surveys — cosmology: observations — X-rays: 
j_j ■ galaxies — X-rays: general. 

1. Introduction 

X-ray emission in galaxies has at a basic level two origins: gravitational potential energy 
release from accretion processes and thermal emission from hot gas. Accretion processes may be 
divided among three main systems: where a lower mass (<1 Mq), and correspondingly longer- 
lived (> 1 Gyr), star accretes onto either a neutron star or black hole (low-mass X-ray binary; 
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hereafter LMXB), where a higher mass (>5 M Q ) and thus shorter-lived (< 10 7 yr) star accretes 
onto a neutron star or black hole (high-mass X-ray binary; hereafter HMXB), and finally where 
accretion occurs onto a nuclear supermassive (10 6 -10 8 M©) black hole (active galactic nucleus; 
hereafter AGN). Hot gas is produced through a combination of stellar mass loss and supernovae, 
and this phase is often the dominant one in the interstellar media (ISM) of galaxies. Indeed, one of 
the more notable results in X-ray studies of galaxies was finding the expected interstellar-medium 
(ISM) component of elliptical galaxies (Forman et al. 1979). The ISM in ellipticals is largely invisible 
at radio and optical wavelengths and had for a time presented a mass-conservation problem given 
the expected output from stellar winds. The hot ISM in galaxies with younger stellar populations 
may be highly extended, as supernova and massive stellar winds give rise to large-scale outflows 
(superwinds; Heckman et al. 1990). For the moment ignoring AGN, it is thus expected that in 
early-type galaxies, with older stellar populations, we will find that X-ray emission is dominated 
by LMXBs and hot ISM. In late-type galaxies, we expect a mix of emission from hot gas, LMXBs, 
and HMXBs with starburst galaxies having a great amount of both hot gas and HMXBs. The hot 
gas component is typically X-ray soft (dominating at < 1 keV) whereas the X-ray binary emission 
is more X-ray hard (dominating the total X-ray emission at > 2 keV). 

It appears that, as expected, the properties of X-ray binaries are correlated with star formation 
rate (hereafter SFR) in many nearby ( < 50 Mpc) galaxies studied intensively with Chandra (e.g., the 
X-ray Luminosity Functions of the binaries appear to obey scaling laws based upon SFR; Kilgard 
et al. 2002; Colbert et al. 2004). The total X-ray emission from binary systems may be written 
as a linear combination of the young population associated with SFR and the older population 
associated with stellar mass (i.e.,Lx = a * SFR + (3 * M; David et al. 1992; Colbert et al. 2004). 
In galaxies with a higher ratio of SFR to stellar mass, where the LMXB component is expected to 
be negligible, the total X-ray luminosity may be dominated by HMXBs and it may be possible to 
use the total X-ray luminosity, neglecting other contributions, as an SFR indicator (e.g., Gilfanov 
et al. 2004a). Indeed, there have been indications that the total X-ray luminosity of a galaxy may 
serve as a star-formation diagnostic at higher redshift (z^0.3, <i;J1700 Mpc) (e.g., Grimm et al. 
2003; Bauer et al. 2002). Of cosmological interest, it also appears that X-ray emission from distant 
Lyman Break galaxies may be similar to galaxies in the nearby Universe (e.g., Brandt et al. 2001; 
Seibert et al. 2002; Nandra et al. 2002; Reddy & Steidel 2004). 

The X-ray band holds much promise as an independent probe of SFR. For instance, it has 
been shown that one of the chief uncertainties in the determination of the power-law index of the 
Initial Mass Function (hereafter IMF) arises mostly from unresolved binary star systems (Kroupa 
2001). Theoretical models of the production of accreting X-ray binary systems, which are probed 
through X-ray observations of star- forming galaxies, are sensitively dependent upon properties of 
the binaries in galaxies (e.g. Sipior et al. 2004). It is thus possible that as both the models and 
the empirical X-ray /SFR relations improve we may develop a probe that is particularly sensitive to 
the binary populations in galaxies. Interesting conclusions regarding the evolution of the "binary 
IMF" have been discussed in Norman et al. (2004). 
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Constraining the X-ray/SFR relation requires understanding the X-ray components of galaxies 
that are not directly associated with current star formation. There have been fairly different 
approaches to this problem. For instance, Grimm et al. (2003) select galaxies with high ratios of 
SFR to stellar mass (hereafter SFR/M*) claiming that they can then reject the LMXB component 
more easily. Grimm et al. (2003) claim that the Lx-SFR relation is linear only at higher SFR 
( > 5 M Q yr -1 ) where high-mass X-ray binaries (HMXBs) dominate the X-ray emission in galaxies. 
At lower SFR, they claim that the HMXB Lx-SFR relation becomes non-linear due to small number 
statistics (Gilfanov et al. 2004b). Four other groups have found that the Lx-SFR relation is linear 
from low to high SFR (Bauer et al. 2002; Ranalli et al. 2003; Colbert et al. 2004; Persic et al. 2004). 
Both Persic et al. (2004) and Colbert et al. (2004) include galaxies with lower SFR/M* values and 
make direct accommodation for the LMXB component. Additionally, Colbert et al. (2004) uses 
only the resolved point source population so that the hot gas and nuclear (AGN) component may 
be rejected. Persic et al. (2004) use X-ray spectral fits to reject emission from hot gas and AGN. 
Bauer et al. (2002), whose work is at higher redshift, select optical emission-line galaxies which do 
not show signs for AGN activity. The result is that there is also quite a spread in the normalization 
of the Lx-SFR relation; published values differ by approximately a factor of five. This difference 
appears to lie with the way in which one treats the LMXB and AGN X-ray "contamination" in the 
total X-ray luminosity of galaxies. 

Understanding this relationship is also quite important for the study of distant galaxies, where 
AGN contamination may be much more difficult to constrain due to the dilution of AGN optical 
spectral features by the host galaxy light (e.g., Moran et al. 2002). High-z classifications are also 
more difficult due to the inaccessibility of the Ha emission line, a useful discriminator between 
true absorption-line dominated (passive) galaxies and those which contain highly obscured AGN 
(e.g., Goto et al. 2003). There is even an extreme class of objects referred to as "X-ray Bright, 
Optically-Normal Galaxies" (XBONGS) which have optical absorption-line dominated spectra but 
very high hard X-ray luminosities (Comastri et al. 2002). Any census of AGN activity in the optical 
waveband must be concerned with the frequency of XBONG-type objects. 

Systematic methods for separating X-ray-detected galaxies dominated by AGN emission from 
those dominated by starburst emission have recently seen a breakthrough with e.g., the Bayesian 
analysis of Norman et al. (2004). This work has allowed the first normal/starburst galaxy X-ray 
luminosity function to be constructed at z ~ 0.3 and z ~ 0.7 (Norman et al. 2004). There is clear 
(and expected) evolution in this XLF, which has a lognormal shape, similar to FIR luminosity 
functions. The evolution appears to be consistent with the expected (1 + z) 2 - 7 evolution of SFR. 
As such high-z normal galaxy XLF studies continue to improve, we may hope to detect subtle 
evolutionary differences between e.g., the X-ray and FIR in order to constrain e.g., X-ray binary 
evolution over long timescales. Such work critically depends upon a thorough understanding of 
the relative contributions of X-ray emission associated with SFR, AGN or quiescent populations in 
galaxies. 

There remains a problem: there still is no suitable normal galaxy sample constructed in the 
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nearby (z<0.1) Universe with which to make comparisons with high-redshift X-ray work (e.g., 
Hornschemeier et al. 2003; Norman et al. 2004). The observational expense of studying extremely 
local (d < 100 Mpc) galaxies to sufficient depth to probe both diffuse X-ray gas and the binary 
LF (down to ~ 10 38 erg s _1 , thus reaching neutron-star binaries) is great; exposures of 50 ks are 
typically required. Currently there are no sufficiently large (> 30 galaxies) samples complete in 
any particular observational property (optical magnitude and/or morphology, starburst properties 
like the ratio of SFR/mass, etc.). Also, while the range of SFR sampled in X-ray-SFR studies to 
date has been quite large, it has in some cases been heavily biased towards higher star-formation 
rates (SFR). The studies aforementioned have typical average SFRs of ~ 10 M yr _1 , going up to 
~ 100 Mq yr _1 . This makes probing the possible transition to non-linear Lx-SFR quite difficult. 
Also, there have been heterogeneous methods used to measure stellar mass and SFR. For instance, 
Grimm et al. (2003) were forced to use approximately 10 separate FIR studies of nearby galaxies 
in order to derive SFRs, and for the same galaxies the FIR flux estimate often differed by factors 
of several. Reliable UV, Ha, and radio SFR estimates are not always available for the same galaxy. 
This is expected to change dramatically with well-observed nearby-Universe samples such as the 
Spitzer Infrared Nearby Galaxies Survey (Kennicutt et al. 2003). Finally, there is some concern 
about "distance-distance" biases introduced by the large range in distances sampled by some of 
these studies. For instance, many studies place high-z galaxies (z > 0.8) from the Chandra Deep 
Fields on the same correlation fits as galaxies in the nearby Universe (d < 50 Mpc). However, 
strong bias is only expected when the range in flux is quite small whereas the distance range is 
large; for most studies this is not a serious problem. 

It is thus critical to observe galaxies with moderate SFR values (~ 1-10 M yr _1 ) as well as 
galaxies with no evidence for current star formation. Since the observational challenge at higher 
redshift is separating the type-2 AGN from the normal/star- forming galaxies, low-to-moderate 
luminosity type 2 AGN should also be included. The sample should have uniformly measured 
SFRs and stellar masses and some highly accurate diagnostic of accretion activity. With these 
requirements in mind, we have turned to the Sloan Digital Sky Survey (SDSS) Data Release 2 
(DR2) catalog and the Chandra X-ray Observatory archive to search for a sample of X-ray-detected 
galaxies near z ps 0.1. DR2 is the second major survey-quality data release of the SDSS (April 
2004; Abazajian et al. 2003). 

At z^0.05 the space density of galaxies with i?<18 reaches ~ 100 deg -2 , so that one would 
expect ss 7 galaxies in each Chandra ACIS-I field of view (« 16' x 16'). The SDSS is covering 
pa 10,000 square degrees to this optical depth. The typical SDSS galaxy redshifts (z ps 0.1) are local 
enough to provide a calibration sample for the high-redshift studies of galaxies. We are working 
towards providing the missing calibration of the X-ray properties of galaxies at z ~ 0.1 by studying 
a uniform sample of SDSS galaxies with serendipitous Chandra detections. 

The major advantage of the SDSS is its uniformity; all galaxies have been observed with the 
same telescope and spectrograph; there is no ambiguity about optical filters or other calibration 
issues. The spectroscopic coverage of galaxies in the SDSS DR2 over its ps 2,600 square degree 
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area is complete down to R ss 17.8 (Stoughton et al. 2002). Chandra has the sensitivity to reach 
galaxies with log {&) ~ —2 (typical of e.g., luminous dusty starbursts, Alexander et al. 2002) in 

J R 

the 0.5-2 keV band at r ~ 18 with 40 ks of ACIS exposure. 1 Thus the depths of coverage of many 
Chandra ACIS observations should be sufficient to place significant constraints on the galaxies 
with SDSS spectra. 

Our work on individually X-ray detected SDSS galaxies is quite complementary to existing 
work. There have been some initial promising results on spiral galaxies from statistical analysis 
of galaxies (not individually detected) using a combination of shallow ( < 5 ks) XMM — Newton 
observations and the 2dF redshift survey (Georgakakis et al. 2003). Georgakakis et al. (2004) 
have also carried out a focused program to identify optically bright, X-ray faint objects which 
are typically quiescent galaxies (e.g. Hornschemeier et al. 2003), reporting on constraints on X- 
ray galaxy number counts at bright X-ray fluxes using SDSS and XMM. Our study only selects 
against Type-I AGN and unresolved sources, we make no exclusion of sources with higher X-ray- 
to-optical flux ratios (which may include some dusty starburst galaxies and higher-mass elliptical 
galaxies). Our focus on optically selected, individually X-ray detected galaxies also complements 
other serendipitous Chandra programs which are driven by X-ray selection and naturally favor 
AGN (e.g., the Chandra Multiwavelength Project, ChaMP; Green et al. 2004). 

The purpose of this paper is to compare optical activity diagnostics such as emission-line 
ratios with the X-ray properties of galaxies. The ultimate goal of our studies will be to assemble 
a complete field galaxy sample; this paper is a pilot study which works towards that ultimate 
goal (in fact, although we avoid cluster centers, a significant fraction of the galaxies are in cluster 
environments). Our chief emphasis is on the X-ray/SFR relation, but we also explore the X-ray 
properties of quiescent galaxies and AGN. A future paper will cover the normal galaxy SDSS XLF. 

The cosmology assumed in this paper is Hq = 70 km s~ l Mpc -1 , Qm = 0.30, S7a = 0.70. 

2. Sample Construction 
2.1. X-ray data 

Our primary goal is to characterize the X-ray emission from star-forming (non-AGN) galaxies in 
the SDSS DR2 database. To accomplish this, we require Chandra sensitivity to detect such galaxies 
at the median SDSS redshift of z ~ 0.1. We only consider Chandra observations of at least 30 ks, 
which ensures that we reach 0.3-8 keV X-ray luminosities of « 2.7 x 10 40 erg s -1 (assuming T = 1.8). 
Based on the Ranalli et al. (2003) SFR/X-ray relation, which is for the total X-ray luminosity, this 
would correspond to a SFR sensitivity of at least 6 M Q yr _1 . Chandra observations significantly 
shorter than this would not be useful for our purposes. 



x The r magnitude here is the SDSS r-band filter. 
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We queried the Chandra archive for ACIS non-gratings observations of at least 30 ks depth 
as of May 2004. We then used various SDSS cross-correlation and finding chart tools (publicly 
available through the main SDSS website 2 to ascertain that 53 fields fell within the SDSS DR2 
spectroscopic footprint. Figure 1 gives the distribution of exposure times for these 53 fields; 24 of 
the fields had X-ray detected SDSS galaxies in them. 

These 53 X-ray datasets were reduced using the XASSIST 3 suite of software tools (Ptak & 
Griffiths 2003). XASSIST uses the basic data reduction steps recommended by the CXC "threads." 
The data are reprocessed to level 2 event files starting from the level 1 events and useful auxiliary 
products such as exposure maps are created. The optional 0.5 pixel position randomization was 
not performed. After reprocessing, each CCD is treated as a different detector and individual CCD 
data are processed separately. For data processing, PI values were constrained to be in the range 
14-548, corresponding to the energy range 0.3-8.0 keV. 

Sources are detected using the CIAO routine WAVDETECT (Freeman et al. 2002), with the 
detection probability set at 1 x 10~ 6 . Background flares were removed by fitting the background 
light curve with a constant and excluding times where the background level deviated by more than 
10<7 from the mean (such a conservative limit was chosen since the background is typically negligible 
for point-like Chandra sources). An image extracted for each source was then fit with an elliptical 
Gaussian model plus sloping background model using Poissonian statistics (i.e., the Cash "C" 
statistic). This provides an improved estimate of the extent and total count rate. The significance 
of each source is then re-assessed based on the total source counts and local background derived 
from the spatial fit. Sources with Poisson probabilities below the 95% level are deleted from the 
source list. Response matrices are then created for each source to estimate the source flux assuming 
a power-law model with T = 1.8 and the Galactic Nh value at the local of the source (determined 
using the FTOOL "nh"). We then visually inspected the sources to look for problems such as 
proximity to CCD chip gaps (resulting in some re-processing) and/or proximity to bright X-ray 
sources (a few sources were rejected). 

As a sanity check, we compared the X-ray catalog produced by XASSIST with the published 
X-ray catalog for the Chandra Deep Field-North (CDF-N Alexander et al. 2003), and found that 
of the 160 X-ray sources identified by XASSIST, 138 have matches in the deeper CDF-N catalog. 
The ones that do not have matches in many cases show signs of X-ray variability with higher 
S/N in the individual shorter observations. The astrometry in the two catalogs matched to within 
~ l'/5. This offset may easily be explained as the difference between the astrometric corrections of 
Alexander et al. (2003) and our "raw" Chandra coordinates. 



2 The SDSS SkyServer website is http://cas.sdss.org/. 
3 See http://xassist.pha.jhu.edu 
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2.2. Optically-derived Galaxy Parameters from SDSS DR2 

The X-ray observed galaxies in this sample were drawn from the SDSS DR2 'main' galaxy sam- 
ple (Strauss et al. 2002). Details concerning the spectral characterization (emission-line measure- 
ment, etc.), stellar mass estimation, and star-formation rate estimates may be found in Kauffmann 
et al. (2003b), Kauffmann et al. (2003a), Brinchmann et al. (2004b) and Tremonti et al. (2004). 
These high-quality measurements were made for 211,894 galaxies and are now publicly available 
(a synopsis of the data products and of the project is given in Brinchmann et al. 2004a, hereafter 
"the SDSS DR2 sample"). This group made made a number of improvements over the standard 
SDSS spectroscopic pipeline in e.g., emission-line flux measurements from continuum-subtracted 
spectra. 

The standard emission-line ratios (e.g., [Nil] /Ha and [OIII]/H(3) were used to characterize 
the galaxy spectra; such line ratio diagrams are often referred to as "BPT diagrams" after Baldwin, 
Phillips, & Terlevich (1981). Higher values of [Nil]/ Ha (> 0.6) generally indicate AGN activity, 
with [Oil I]/ 'HP discriminating between Low Ionization Nuclear Emission- line Regions (LINER; 
Heckman 1980) and Seyfert-type objects (LINERs have [Oil I]/ HP < 3). Kauffmann et al. (2003a) 
have suggested two lines of demarcation (see Figure 2) between AGN and starburst galaxies in the 
BPT diagram which broadly separate galaxies into three classes: AGN, star-forming galaxies, and 
AGN/star- forming composites. The top line is the theoretical upper bound for starburst galaxy 
emission-line ratios from Kewley et al. (2001). Galaxies with values to the right (above) this line in 
the BPT diagram are expected to have appreciable contribution from an AGN to their emission-line 
flux. The second boundary, below which AGN activity is not expected and galaxies are typically 
star-forming galaxies, is a modified version of the Kewley et al. (2001) line given by (Kauffmann et 
al. 2003a; see Figure 2). The optical classification of galaxies between these lines is AGN-starburst 
"composite". Composite galaxies are expected to have ss 41% of their [OIII] luminosity density 
and ss 11% of their Ha luminosity density arising from an AGN (Brinchmann et al. 2004b). 

One problem encountered in any star formation study is that the SFR cannot be measured 
directly. A large number of star formation indicators have been theoretically or empirically cal- 
ibrated at wavelengths ranging from the UV to the radio. In this work we adopt the total (i.e. 
aperture-corrected) SFRs of Brinchmann et al. (2004b). These SFRs are derived from the optical 
nebular lines (principally Ha) for galaxies classified as star-forming, and from the 4000A break for 
galaxies classified as AGN. The majority of previous X-ray/SFR studies have relied on far-infrared 
(FIR) SFR estimates. However, Kewley et al. (2002) demonstrated that Ha and FIR SFR estimates 
agree to 10% when proper account is made for attenuation. 

For a full description of the SFRs the reader is referred to Brinchmann et al. (2004b); we 
briefly summarize the methodology here. The star formation rate within the SDSS fiber apertures is 
estimated from the strong optical emission lines ([Oil], [OIII], H/3, Ha, [Nil], [SII]) using a Bayesian 
approach and a grid of 2 x 10 5 nebular models (Chariot & Longhetti 2001). The SFR is estimated 
simultaneously with the metallicity, ionization parameter, dust attenuation, and dust-to-metal ratio 
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in order to account for variations in the unattenuated Ha/H/3 ratio and the conversion between Ha 
luminosity and SFR. This approach differs from the standard Kennicutt (1998) translation between 
Ha and SFR but there is good agreement in the mean conversion factor (see Figures 7 and 8 of 
Brinchmann et al. 2004b). For galaxies harboring an AGN, the SFR in the fiber is estimated from 
the D4000 spectral index, a measure of the 4000 A break. This secondary SF indicator is calibrated 
consistently from the star-forming galaxy sample. Because the SDSS fibers only encompass ~ 1/3 
of the galaxy light, the SFRs must be aperture corrected. The galaxy colors outside of the fiber 
aperture are used to estimate this correction. Brinchmann et al. (2004b) adopt a Kroupa IMF 
(Kroupa 2001). To convert to a Salpeter (1955) IMF between 0.1 and 100 M Q , one would multiply 
the Brinchmann SFRs by a factor of 1.5. In this paper we correct the SFRs of our comparison 
sample (see §2.3) to the Kroupa IMF. 



2.3. SFR measurements of the comparison studies 

There are five studies which we compare in this work (in order of discussion): Colbert et al. 
(2004), Persic et al. (2004), Bauer et al. (2002), Ranalli et al. (2003) and Grimm et al. (2003). All 
of the SFR measurements may be traced back to the FIR-SFR calibration of Kennicutt (1998), 
which is based upon the "standard" Salpeter IMF discussed above. However, in this section we 
aim to quantify the differences and uncertainties in order to uncover any real disagreement among 
the various X-ray/SFR relations. 

Colbert et al. (2004) and Persic et al. (2004) used Kennicutt (1998) FIR-SFR conversions. 
However, both made some additional correction to the IRAS FIR luminosity. Colbert et al. (2004) 
added the ultraviolet luminosity to the IRAS FIR luminosity to produce a "corrected" FIR lumi- 
nosity (which is higher than the directly measured FIR by up to a factor of 2) and then apply a 
bolometric correction of a factor of 1.4 (thus their FIR luminosities may be up to 40% higher than 
IRAS FIR luminosities). Persic et al. (2004) note that at lower SFR, non-SF-related emission such 
as cirrus may be present and use the i?-band luminosity to make this correction (a "minor" 14% 
correction). Therefore, both these studies (Colbert and Persic) are based upon IRAS fluxes with a 
14-40% correction factor applied. 

Ranalli et al. (2003) used the Kennicutt (1998) FIR-SFR conversion (based upon the standard 
Salpeter IMF) and the Condon (1992) radio-SFR conversion. Inspection of the Ranalli equations 
indicates consistency with the Kennicutt (1998) calibration despite the Condon (1992) radio-SFR 
conversion factor being for high- mass stars (> 5 Mq) rather than for the more inclusive mass range 
for the Salpeter IMF (0.1-100 M ). 

Bauer et al. (2002) use the radio-SFR relation of Yun et al. (2001), who derive the FIR-radio 
correlation and then use the Kennicutt (1998) FIR-SFR relation to derive a radio-SFR conversion. 
The Yun et al. (2001) SFR-radio normalization is lower than the Condon (1992) normalization in 
the sense that Yun predicts less radio emission per unit SFR by a factor of 2.3. Since the Yun 
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SFRs are for the standard Salpeter IMF over a 0.1-100 M interval but the Condon (1992) SFRs 
are for a narrower mass interval (> 5 MQonly), this difference in the calibration of the 1.4 GHz 
estimator is as expected. Thus, Bauer et al. (2002) uses a 1.4 GHz estimator that appears consistent 
with the rest of the FIR measures and with the standard Salpeter IMF. 

Grimm et al. (2003) is somewhat more difficult to compare with the others. The FIR-SFR 
conversion used is the same Kennicutt (1998) conversion used by the others, but their measurements 
are much more heterogeneous than the other authors. Grimm et al. (2003) have searched the 
literature for SFR measurements in nearby galaxies at Ha, ultraviolet, FIR, and radio wavelengths. 
The Ha-SFR conversion they use, from Rosa-Gonzalez et al. (2002), is 40% higher than Kennicutt 
(1998) in the sense that the derived Grimm SFRs would be higher using the Rosa-Gonzalez et al. 
(2002) value. The radio conversion they use is from Condon (1992), but they use the thermal radio- 
SFR conversion factor rather than the non-thermal (which are different by an order of magnitude). 
There appear to be relatively few radio measurements in the Grimm paper, so this appears not to 
matter. Thus it appears that Ha measurements could result in SFR measurements that are up to 
40% higher than the other authors, but these are not always used in their estimates. 

2.4. Chandra-SDSS source matching 

In the 53 ACIS fields there are 338 "well-observed" SDSS DR2 galaxies with 0.03 < z < 0.25. 
Our definition of well-observed is that Chandra covered at least 95% of the galaxy's 50%-encircled 
light radius (we use the Petrosian 50% encircled-light radius in the g-band) and that the galaxy 
was observed at a Chandra off- axis angle < 12'. Including galaxies with off-axis angles larger than 
12' resulted in few additional detections. 

The SDSS astrometry is excellent (< Of! 5) and as mentioned in § 2.1, the X-ray positions from 
XASSIST are accurate to better than 2"0. The optical extent of the SDSS galaxies was typically 
larger than either of these two positional errors (the median galaxy radius was « 2 / .'2; see Figure 3a). 
SDSS galaxies were considered matched to an X-ray source if the X-ray source position was within 
the optical 50% encircled-light radius (we used the Petrosian radius in the g-band). We found 
that 42 galaxies contained X-ray detections and the X-ray detections were typically well within the 
optical galaxy extent (see Figure 3a). None of these 42 sources were the primary target of the X- 
ray observation but there is clearly some bias towards cluster members, as indicated in the redshift 
distribution in Figure 3b. We randomly offset the SDSS galaxy source positions by 20" and 30" 
to evaluate the false-match probability and found that < 1 false match is expected among the 42 
galaxies. The X-ray and optical properties of these 42 galaxies are given in Table 1, including the 
X-ray fluxes, redshifts, absolute magnitudes, and optical spectroscopic classifications. The X-ray 
source SDSS spectra and images are shown in Figure 4. The median number of 0.5-8 keV X-ray 
counts per source is just 18.2, so we do not present X-ray spectroscopic results in this paper. 

We have examined these sources carefully to ascertain that they are not affected by instru- 
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mental problems such as proximity to a chip gap in the ACIS camera and/or proximity to any 
bright X-ray source. The Galactic column density along the various lines of sight to these galax- 
ies are fairly low, the maximum Ah is = 4.7 x 10 20 cm" 2 , the mean and standard deviation is 
Ah= 2.6 ± 1.2 x 10 20 cm" 2 . X-ray fluxes are corrected for Galactic absorption. 

There are published results on 37 X-ray-detected normal and starburst galaxies in the CDF-N 
(Hornschemeier et al. 2003) so we again used this field to test our results. These 37 CDF-N X-ray 
detected galaxies are typically at higher redshifts and thus below the optical flux limit of SDSS 
(median R 19 mag). Three of the four SDSS DR2 galaxies in this field were X-ray undetected 
in our sample (which did not use the full combined depth of the CDF-N). One additional source 
is detected in the full 2 Ms depth of (Hornschemeier et al. 2003) and the final two were at large 
off-axis angles where the effective exposure time is low. This illustrates that the extreme depth 
of the CDF-N is not necessary for this kind of lower redshift work and also shows the necessity of 
building up large areal coverage to probe the relatively sparse galaxy population at z ~ 0.1. 

3. Results 

3.1. Optical Galaxy Classification 

The majority of the Chandra-detected SDSS DR2 galaxies are "emission-line" galaxies but 
the precise number depends upon the significance cutoff used for the emission line flux. Note that 
all classifications are based upon continuum-subtracted spectra and that the optical spectroscopic 
apertures typically cover ~ 30% of the host galaxy. 

The emission-line ratios for both the X-ray detected galaxies and the full SDSS DR2 galaxy 
sample are shown in Figure 5a and their classifications based upon those ratios are given in Table 2. 
Galaxies are classified if they have all four emission lines ([OIII], H/?, [Nil], Ha) detected at > 2a; 
there are 24 galaxies matching this criteria (see Table 2). It is also possible to use the [Nil]/ Ha 
ratio alone (see Figure 5b) as a coarse AGN discriminator; this yields an additional 3 emission-line 
galaxies, all of which are classified as AGN. Note that at high-redshift, Ha would be unobservable 
so these galaxies would likely have been misclassified. 

There are thus 24 emission-line galaxies (ELGs), three lower-significance ELGs (with only [Nil] 
and Ha detected above 2a), and 15 galaxies which appear to be absorption-dominated to fairly 
strict limits. 

3.2. X-ray/SFR Correlations 

In Figure 6 we show the total integrated X-ray luminosity in the 0.3-8 keV band for the 
Chandra-detected ELGs plotted against SFR as measured by Brinchmann et al. (2004b). We also 
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show the correlations between SFR and X-ray emission found by several other authors. We have 
adjusted these various correlations to the 0.3-8 keV bandpass. We assumed a T = 1.8 power-law 
for conversion; the correction is a factor of ~ 2 from the 2-8 keV bandpass to the 0.3-8 keV 
bandpass. In Section 2.3 we describe the various SFR measurements used by these other authors; 
their methods yield differences of up to 50% in SFR. In Figure 6 we divide the SFRs from the 
comparison studies (Colbert, Persic, Grimm, etc.) by a factor of 1.5 to account for the Kroupa 
IMF. The Colbert SFRs were additionally divided by a factor of 1.4 to remove the bolometric 
correction over SFRs derived purely from IRAS fluxes. 

There are seven pure "star-forming" galaxies in our sample based upon emission-line ratio 
diagnostics alone. There is one clear outlier among the X-ray-detected "star-forming" galaxies in 
Figure 6 (X102348. 48+040552. 4; it is marked in cyan in the BPT diagram of Figure 5) that has 
a much higher X-ray luminosity than is expected based upon its SFR. Inspection of its optical 
spectrum reveals that it is a Narrow-Line Seyfert 1 (NLS1; see Figure 4). If one compares the 
width of the Balmer lines in this spectrum with the widths of the forbidden lines, the Balmer width 
is found to be much larger. Thus while optical emission line ratio diagnostics are fairly good at 
separating AGN from star-forming galaxies, we find that the additional information afforded by line 
widths is quite critical. We have plotted the Balmer vs. forbidden line widths for all the emission- 
line galaxies in Figure 7. As can be seen, there are only three other sources with Balmer widths 
that are appreciably larger than the forbidden line widths; these three are composite galaxies (or 
deemed to be low signal-to-noise AGN). 

Correcting for the likely LMXB contribution using the stellar mass estimates from Kauffmann 
et al. (2003a) and the LMXB mass constant of Colbert et al. (2004) brings five of the six remaining 
star-forming galaxies into agreement with the Persic et al. (2004) line; one of these has such a 
low SFR that all of its X-ray emission is explained by LMXB emission (it thus does not appear 
in the second panel of Figure 6. There remains only one "star-forming" galaxy besides the NLS1 
that has a higher value of X-ray emission than expected from the Colbert relation. This galaxy, 
X150944. 45+570434. 5, is in a particularly clean Chandra field (there is a z ss 4 QSO in the field, 
but no large scale structure, etc.). The source has no detected 2-8 keV counts and thus does not 
show evidence for harboring an AGN. 

3.3. Galaxies with Absorption-line Dominated Optical Spectra 

There are 15 galaxies which lack emission lines at the la level. Expecting that these galaxies 
do not have much current star-formation, we plot their total X-ray luminosities versus their host 
galaxies' stellar masses (Figure 8; the masses are the medians of the probability distributions of 
dust-corrected mass for each galaxy from Kauffmann et al. 2003b). The total X-ray luminosities 
exceed that expected from X-ray binaries alone. The Lx-mass relation is consistent with the Lx~ 
Lb relation of Fabbiano et al. (1992) which was attributed to the X-ray emission from hot gas in 
galaxies. 
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The spectrum of the one outlier in the absorption-dominated galaxy X-ray/mass plot, X171206.83 
+640830.7, does not show any sign of AGN emission (see Figure 4) and has been previously re- 
ported as a post-starburst galaxy by Davis et al. (2003). The X-ray emission is extremely soft 
(r ss 3.4, see Figure 9) and the X-ray luminosity is quite high (Lx ~ 1-6 x 10 42 erg s _1 ) but the 
source is detected in the 2-8 keV band as well with Lx ~ 1.3 x 10 41 erg s _1 (2-8 keV). There is 
no detected SFR within the aperture (no Ha emission) and the estimated current SFR for the 
entire galaxy is 0.6 M yr -1 , with a 97.5% upper limit at 4.1 M yr _1 (Brinchmann et al. 2004b). 
The contribution from either an evolved component or current star-formation to the hard X-ray 
luminosity is estimated to be < 5%. Davis et al. (2003) indicate that this source may be variable, 
having declined in soft X-ray flux by a factor of several since ROSAT observations in 1993 and 
1994; they therefore suggest that this source is either an AGN or a bright Ultra-Luminous X-ray 
source (ULX; e.g., Colbert & Ptak 2002). Our continuum-subtracted optical spectra place strong 
constraints on the AGN nature of this object (we do not detect any AGN signature) and if the 
lack of emission lines were due to obscuration, we might have expected to see evidence in the X-ray 
spectrum, which we do not (see Figure 9). So we explore the post-starburst nature of this object 
in more detail. 

Kauffmann et al. (2003b) calculate burst mass fractions for such galaxies. The median of the 
mass fraction distribution indicates that 10% of the mass was produced in a burst, corresponding 
to f» 2 x 10 9 M Q . However, as noted in Kauffmann et al. (2003b), there is an age/mass degeneracy: 
a large burst that occurred long ago is indistinguishable from a smaller burst that occurred more 
recently. The distribution of the mass fraction actually extends from zero to 65%, the zero case 
corresponds to the stellar mass having been formed a long time ago (> 2Gyr). We adopt the 
median value of 10% acknowledging that there is some significant uncertainty. 

Sipior et al. (2004) model the X-ray binary population resulting from a 10 M yr _1 starburst 
lasting 20 Myr; the 2-10 keV luminosity of such a burst has a peak near 2 x 10 40 erg s- 1 10 Myr 
after the burst has ended (assuming a Salpeter IMF), declining to < 10 38 erg s _1 by 1 Gyr after 
the burst. If X171206. 83+640830. 7 were a post-starburst galaxy with 10% of its mass formed in 
a 20 Myr interval, the implied star-formation would be 150 M Q yr _1 (scaling by a factor of 1.5 
to account for the Kroupa IMF). Assuming we observe X171206. 83+640830. 7 very soon after the 
burst (near the X-ray maximum at burst+10 Myr) then scaling the Sipior et al. (2004) curves yields 
an SFR of ~ 65 Mq yr _1 , which within errors is very consistent. We find the coarse agreement 
between likely post-starburst nature of this galaxy and the X-ray emission encouraging. 

However, the lack of emission lines indicate that this source must be viewed much more than 
10 Myr after its peak in X-ray luminosity. Since the peak X-ray luminosity would be even higher 
than observed now, the implied evolved starburst strength would be scaled upward by nearly an 
order of magnitude. A more plausible explanation for the high X-ray luminosity, given the X-ray 
variability, is that this is an extremely X-ray luminous accreting binary. Such a source would thus 
be the most X-ray luminous Ultraluminous X-ray source (ULX; see e.g., Fabbiano et al. 2003) ever 
detected. 
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3.4. AGN Activity and XBONGS 

There are 19 galaxies in this sample which may be classified as AGN (see Table 2). Their 
Lx/£[om] ratios span the range found for nearby Seyfert 2 galaxies (e.g., Pappa et al. 2001). 
About 12 of the galaxies are at the upper end of the expected range, indicating possible Compton- 
thick obscuration (Lx/L\oiir\ > 10). The statistics on the hardness ratios for these 12 objects 
[ratio of (H-S)/(H+S) where H=2-8 keV counts and S=0.5-2 keV soft counts] are not sufficient 
to test for a correlations with Lx/^ro/jn ( as was done for a nearby sample of Seyferts in Pappa 
et al. 2001). We will be exploring these high Lx/L^qiii] AGN in more detail in a future paper (A. 
Hornschemeier et al., in preparation). 

There are three AGN with no detected [OIII]/H/3 but Ha/ [Nil] ratios indicating they are 
AGN. These three sources are candidate "XBONGS" (X-ray Bright, Optically Normal Galaxies) 
which is an observational class of AGN found in deep X-ray surveys with no optical signature 
for AGN activity. Similar objects at higher redshift (in deeper surveys) would likely have been 
classified as absorption-dominated galaxies due to the inaccessibility of the Ha and [Nil] emission 
lines. These three candidate XBONG galaxies do not, however, have high X-ray luminosities 
(L x < 2 x 10 41 erg s _1 , 0.3-8 keV). Thus they do not appear to be the observational analogs to 
XBONGS found in hard X-ray surveys (e.g., Comastri et al. 2002). There are similarly no extremely 
hard, X-ray luminous sources among the 15 X-ray detected absorption-dominated X-ray-detected 
galaxies. There are a total of 103 absorption-dominated galaxies in all of the Chandra fields at 
off-axis angles < 8' (13 of the 15 X-ray detected absorption-line galaxies are at off-axis angle < 8'. 
Thus we do not find any galaxies which show signs for extremely luminous AGN in the X-ray band 
that show no optical signature of activity in the high signal-to-noise SDSS spectra. 



4. Conclusions 

We have conducted a pilot study of 42 SDSS DR2 galaxies serendipitously detected in the 
X-ray band. The high signal-to-noise SDSS spectra afford detailed activity diagnostics that are 
often missing in the deepest surveys. 

Among X-ray detected galaxies which show no signs for AGN activity (15 absorption-dominated 
galaxies and 7 star-forming galaxies), the X-ray emission is naturally explained by a combination of 
LMXBs, HMXBs, and hot ISM. Thus, we do not find evidence for AGN that are missed in the SDSS 
spectroscopic analysis. Our work thus supports that of others (e.g., Moran et al. 2002) that fairly 
simple observational effects may make X-ray detected AGN appear as XBONGS at higher-redshift. 
Among the 19 SDSS galaxies which have AGN signatures in their optical spectra, the Lx-£[om] 
ratios are found to be coarsely consistent with Seyfert 2's in the nearby Universe. 

We do find four instances where galaxies might have been misclassified in lower signal-to-noise 
optical spectra. A NLSl's line widths clearly reveal the presence of an AGN while its emission-line 
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ratios are similar to that of a starburst galaxy. Three other galaxies have pure absorption-dominated 
optical spectra blueward of Ha but exhibit Ha and [Nil] emission typical of AGN. All three of these 
galaxies would not have been identified as AGN above z ss 0.5. 

The detection of X-ray emission from a post-starburst galaxy indicated that the burst mass 
fraction calculated from the optical spectroscopic models may be in agreement with the expected 
X-ray emission from a starburst. However, the likely evolved nature of the post-starburst (the age 
of which is poorly constrained) indicates that the X-ray luminosity is quite excessive given the 
magnitude of the starburst. This one object is not sufficient to test the models but is strongly 
suggestive that a large sample of post-starburst galaxies (perhaps with stronger constraints on the 
burst mass fractions) may have a strong impact on the modeling of accreting binary systems. It 
may also be an excellent way to search for young ULX sources as this galaxy appears to harbor a 
variable X-ray source. We are able to rule out that this post-starburst galaxy is an AGN to very 
strong limits from optical spectroscopy. 

We have made an effort to compare the various Lx-SFR relationships in the literature. We 
find that differences among e.g., bolometric corrections and IMF account for up to a 50% difference 
in the SFR relationships. Any difference greater than this level must then be attributed to e.g., 
a different treatment of the contributions of AGN vs. star-formation. Our analysis supports the 
lower normalization of the Lx-SFR relation after proper account of AGN activity and contributions 
from evolved stellar populations are made. 
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Table 1. SDSS DR2 Galaxies Detected Serendipitously by Chandra 
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X011313. 86+152938.2 
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"Star- forming" , NLSl 
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Absorption 
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X163629. 56+410221. 6 


1.4 
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39.0 


0.047 
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10.5 


Composite 
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-23.04 


1.30 


8.04 


11.2 


Composite 


X171206.83+640830.7 


1.1 


402.1+21.0 




96.5 


0.082 


-21.01 


-0.21 




10.3 


Absorption 


X171232.43+640008.7 


1.3 


16.6+4.5 




3.8 


0.085 


-22.82 


-0.51 




11.3 


Absorption 


X171254.28+635930.2 


2.6 


7.2+3.2 




1.7 


0.083 


-23.13 


-0.40 




11.3 


Absorption 


X171303.84+640700.7 


1.1 


23.7+6.2 




6.1 


0.081 


-22.85 


-0.28 




11.4 


Absorption 


X171329. 16+640247.6 


1.1 


71.3+10.4 




17.0 


0.078 


-23.44 


0.13 




11.5 


AGN(< la) 


X221722. 73+002107.0 


1.7 


17.1+4.7 




2.2 


0.095 


-22.12 


1.00 


6.39 


11.0 


Composite 


X234817.93+010615.2 


2.1 


18.7+6.4 




7.8 


0.093 


-23.30 


-0.39 




11.3 


Absorption 



a Absolute magnitude in the 2-band from Kauffmann et al. (2003a) 
b Star- formation rate from Brinchmann et al. (2004b) 

c The median of the dust-corrected stellar mass distribution from Kauffmann et al. (2003b) 

d This source has emission-line ratios typical of a star-forming galaxy but it is in fact a NLS1. See §3.2. 
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Table 2. Optical Spectroscopic Type of Serendipitous Chandra Detections 



Optical spectral type a 


Number of X-ray detected galaxies 




> 3a class b 


> 2a class b 


Ha only c 


Star-forming 


7 


7 




Composite 


8 


9 




LINER 


5 


6 




Seyfert 


1 


1 




Total Emission Line: 


21 


24 


27 c 


Absorption 


21 


18 


15 


TOTAL 


42 


42 


42 



a As determined from emission- line ratios only (see Figure 5. 

b Classifications made for objects who had all four emission- lines 
measured above this significance. Thus, the "absorption-dominated" 
galaxies are those galaxies for which emission lines may be present, 
but at least one is below the Nc value. 

c Emission-line classification based on only the log [Nil]/ Ha ratio. 
13 sources are AGN and 15 sources are star-forming galaxies. 
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20 40 60 80 100 120 140 

Chandra Exposure Times (ks) 



Fig. 1. — X-ray exposure times for the 53 Chandra ACIS fields in this study. There are three fields 
at or below the 30 ks lower cut-off because the net effective exposure time was slightly shorter than 
the nominal observation time (which was used to query the archive). 
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Fig. 2.- 



Emission-line ratios for the total SDSS DR2 sample of Brinchmann et al. (2004a). 
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Fig. 3. — (Left) Comparison of the optical extent of the galaxies with the X-ray/SDSS offsets. 
Both are in units of arcseconds. (Right) Redshift distribution of the 42 serendipitously-detected 
SDSS DR2 galaxies in Chandra ACIS fields (black line) as compared to the full X-ray observed 
sample of 338 SDSS DR2 galaxies (filled gray histogram). The full X-ray observed sample has been 
renormalized for ease of comparison. The properties of these 42 galaxies are given in Table 1. 
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Atlas figure of spectra 
and images available at : 



http : //www.pha . jhu . edu/~annh/SDSSXray 



Fig. 4. — Optical images and spectra for all the SDSS-Chandra matches. The red circle is the 
approximate position of the fiber aperture and the SDSS Right Ascension and Declination is given 
in the upper left corner of each image. 
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Fig. 5. — (Left) Emission-line ratios for the total SDSS DR2 sample (Gray circles) and the X-ray 
detected sources (black symbols). Empty symbols indicate sources with lower-significance emission 
lines (> 2a) whose classifications are thus slightly more uncertain. One slightly lower significance 



source is off this diagram at 



[0.8,0.4]. The breakdown of 



log {[Nil]/ Ha), log ([01 1 1}/ HP) 

source classifications is given in Table 2. The source marked with a cyan symbol is a Narrow- 
Line Seyfert 1 galaxy. The vertical line indicates the [Nil]/ Ha = 0.6 demarcation for AGN 
vs. starbursts and the horizontal line indicates the [01 1 1]/ H (3 = 3 demarcation for Seyferts vs. 
LINERs. (Right) log ([NII]/Ha) for the 28 X-ray detected galaxies with both [Nil] and Ha 
detected above 2a (this is a gain of an additional 3 emission- line galaxies). There are 13 galaxies 
with [Nil]/ Ha > 0.6, versus 10 when using both [OIII]/Hf3 and [Nil]/ Ha. 
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Fig. 6. — (a) Total Chandra full-band X-ray luminosity versus star-formation rate (SFR) as 
calculated in Brinchmann et al. (2004b). The lines show the correlations derived in several studies 
in the literature, all of these studies have been adjusted to the 0.3-8 keV band assuming a T = 1.8 
power law spectrum and SFRs have been corrected to a Kroupa IMF. We have divided the Colbert 
line by a factor of 1.4 to remove their bolometric correction (see § 2.3). We also show the data on 
nearby galaxies from Persic et al. (2004). The dashed box shows the region of the upper limits for 
the 295 X-ray-undetected SDSS galaxies, which are fairly evenly distributed throughout this area 
(plotting them individually confuses the plot), (b) Same as (a) but now the X-ray emission that is 
not associated with current star-formation has been removed from the full-band X-ray luminosity by 
accounting for the LMXB component through stellar mass estimates. Note that one "star-forming" 
galaxy is not shown in panel (b) because its entire X-ray luminosity is explained by the non-SFR 
component (appear near logLX,logSFR = [40,-0.3] in panel (a)) For the SDSS DR2 datapoints 
we have used the mass/LMXB coefficient from Colbert et al. (2004) and subtracted this value from 
the total Lx- The Persic et al. (2004) data points in this graph are now the "corrected" data points 
in Figure 4 of their paper, which has had a correction factor applied for the LMXB component; 
note that the difference is greatest for the lowest SFR. 
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Fig. 7. — fBaimer versus o"Forbidden for the full SDSS DR2 sample having emission lines detected 
above 2a and for the X-ray detected emission-line galaxies (classification symbols the same as in 
previous figures). The NLS1, X102348. 48+040552. 4, is marked in cyan as in Figure 5. 
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Fig. 8. — Full-band X-ray luminosity versus stellar mass calculated in Kauffmann et al. (2003b) 
for the 15 "pure" absorption-dominated galaxies. These are the galaxies which do not demonstrate 
emission lines in any of the four "BPT diagram" lines (see Figure 5) above la in the continuum- 
subtracted spectra. These are very strong limits as these diagnostic spectra cover the inner 30% 
of the host galaxies and thus suffer much less dilution than do higher-redshift samples. The dark 
solid line indicates the stellar-mass/X-ray 'point source relation of Colbert et al. (2004) which is the 
expected LMXB component in these galaxies. The dashed line indicates a fit to the data, excluding 
the outlier X171206. 83+640830. 7 (marked with a diamond). The slope of this line is consistent with 
the non-linear Ly^-L-Q relation found by Fabbiano et al. (1992), with slope ~ 1.8. The non-linearity 
of the relation was attributed to the hot ISM of elliptical galaxies. 
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X-ray spectrum of the outlier X171206. 83+640830. 7 from Figure 8. 



